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BOUNDARY-LAYER TRANSITION ON THE AVERAGE HEAT TRANSFER 
TO A YAWED CYLINDER IN SUPERSONIC FIOW 
By Ivan E. Beckwith and James J. Gallagher 

SUMMARY 


The average heat-transfer to the forward half of a circular cylinder 
was measured in a blowdown jet at a stream Mach number of 4.15* The test 
Reynolds number based on the diameter of the cylinder was varied from 

5 X lO^ to 13 X 10^ and the yaw angle was varied from 0° to 60°. Most 
of the tests were made at a wall-to-stagnation temperature ratio of 
about 0.8. 

The test results showed that, as the yaw angle was increased from 0° 
to 35°> the average heat-transfer coefficient increased by 50 percent at 

a Reynolds number of 5 X 10^, and at a Reynolds number of 13 X 10 ^ the 
heat-transfer coefficient increased by 100 percent. Comparison of the 
measured heat-transfer coefficients with theoretical predictions and with 
other experimental results indicated that these large increases in heat 
transfer, with increasing yaw angle and fixed Reynolds number, were caused 
by transition from laminar to turbulent boundary layer. This transition 
apparently occurred primarily as a result of a dynamic type of instability 
which is associated with the secondary flow in the boundary layer on the 
yawed cylinder. 

Further increase of the yaw angle from 35° to 60° resulted in a 
reduction in the average heat-transfer coefficients of about 35 percent 
from the peak values . Most of this reduction apparently corresponds to 
the decrease in heat transfer with increasing yaw angle that could be 
expected to occur on a cylinder with a completely developed turbulent 
boundary layer. 
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INTRODUCTION 


The use of sharp leading edges on wings and tail surfaces is unde- 
sirable at flight speeds for which aerodynamic heating becomes a serious 
problem because of the high local heat-transfer coefficients and small 
heat capacity of the sharp edge. Both of these difficulties can be 
alleviated by using a blunt leading edge. Blunting the leading edge, 
however, has the disadvantage of increasing the drag. The increased drag 
of the blunted section can be reduced by increasing the sweepback angle 
of the wing. For these reasons, data are needed on the effects of sweep- 
back on the heat transfer to blunt leading edges. 

A yawed cylinder of high aspect ratio conveniently simulates the 
leading edge of a sweptback wing and also allows a simplified treatment 
of the theoretical problem. Thus, for incompressible flow with a Prandtl 
number of 1, inspection of the laminar -boundary- layer equations for the 
infinite yawed cylinder (refs. 1 and 2) shows that the velocity distri- 
bution in planes normal to the cylinder axis and the heat transfer 
depend only on the flow component normal to the cylinder. Solutions for 
this case (refs. 3 and 4) then show that the local Nusselt number, based 
on cylinder chord measured normal to the leading edge, is proportional 
to the square root of a Reynolds number based on the chordwise component 
of the stream velocity and the normal chord. That is to say, for fixed 
stream conditions and fixed normal chord, the heat— transfer coefficient 
is proportional to the square root of the cosine of the yaw angle. Appli- 
cation of these results to a wing of fixed cross-sectional shape, plan- 
form area, and normal chord indicates that, for given stream conditions 
and wall temperature, the total quantity of heat transferred to the wing 
is also proportional to the square root of the cosine of the sweepback 
angle. 

The results of theoretical and experimental investigations made at 
supersonic speeds (refs. 5 8) show that, for laminar boundary layers, 

somewhat larger reductions in heat transfer can be obtained by yawing the 
cylinder than would be expected for the corresponding incompressible-flow 
case. This is largely because of the effect of the detached bow wave on 
the local flow conditions near the cylinder. Apparently, then, as long 
as the boundary layer remains laminar, some reduction in total heat trans- 
fer to a wing and its leading edge is to be expected by increasing the 
sweepback angle. It could also be expected that, as the flight Reynolds 
number is increased, transition from laminar to turbulent flow will 
eventually occur and cause the associated large increases in heat-transfer 
rates . 

A related phenomenon, which is probably not so generally appreciated, 
is the effect of the sweepback or yaw itself on transition. Recent inves- 
tigations (refs. 9 to 13 and papers of limited availability by the British 
investigators P. R. Owen and D. G. Randall) have shown that sweepback or 
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yaw has large adverse effects on transition on a wing or cylinder. Owen 
and Randall have shown that these effects are caused by a dynamic type 
of instability which arises as a result of the secondary or three- 
dimensional flow that occurs within the boundary layer on a sweptback 
wing or cylinder. They also defined and, from experimental data, evalu- 
ated a stability parameter which has been successfully applied to the 
calculation of a lower limit of Reynolds number and yaw angle at which 
this instability first appears for incompressible flow. The results of 
reference 12 indicate that this instability may appear at a much lower 
Reynolds number on a wing in supersonic flow than predicted by the 
Owen-Randall parameter for incompressible flow. 

The purpose of the present report is to present data showing the 
effect of yaw on the heat transfer to cylinders over a higher Reynolds 
number range than that of previous data. The Reynolds number range of 
the present tests, based on free-stream conditions and cylinder diameter, 

is from 5 X 10^ to 13 X 10 5 . This corresponds to flight conditions at 
altitudes of 45,000 to 62,000 feet at a Mach number of 4 with a wing 
leading-edge radius of 1 inch. The data are compared with correlation 
parameters for laminar and turbulent heat transfer and also with laminar- 
heat-transfer data from previous investigations. Estimated values of the 
transitional Reynolds numbers from the present investigation are compared 
with values predicted by the Owen-Randall stability parameter evaluated 
from the data of reference 9 * 


SYMB0IS 


A cross-sectional area of model 

C constant in equation (l4) 

C]_,C 2 constants in equation ( 10 ) 

Cp "normal" pressure coefficient 

c specific heat of model material 

Cp specific heat of air 

D diameter of model 

h heat-transfer coefficient 

k thermal conductivity of air 
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km thermal conductivity of model material 

k',n.' reference conductivity and viscosity evaluated at T' 


L 

length from leading edge or stagnation point along an external 
streamline 

l 

spanwise length of model 

M 

Mach number 

% 

normal component of Mach number ahead of cylinder 

Mi 

local chordwise component of Mach number 

Nu 

Nusselt number 

P 

periphery of model exposed to stream 

* 

P 

pressure 

% 

convective heat-transfer rate per unit area at the surface 

R 

Reynolds number 

r 

recovery factor 

T 

temperature , °R 

T 

e 

equilibrium temperature, °R 

T' 

reference temperature, i^T g + ^w) 

t 

time 

U 

resultant velocity ahead of cylinder 

U N 

normal component of velocity ahead of cylinder 

u 

local chordwise velocity 

V 

local resultant velocity 

V 

local velocity component normal to external streamline • 
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y 

z 

7 

A 

5 

0 


M- 

P 


cr 

T 

X 

♦ 

0) 


distance from stagnation point along the surface of the 
cylinder in a plane normal to the leading edge 


distance measured along span of model 
boundary-layer coordinate normal to surface 
ratio of specific heats 
boundary-layer momentum thickness 
boundary-layer thickness 

angular distance around cylinder measured from stagnation line 

du n 
P 

Pohlhausen boundary-layer parameter, — 

dynamic viscosity of air 
mass density of air 


dx 




T„ 


reference density, — )1 + — )p_ 

2 \ V s 

Prandtl number 
shear stress 
stability parameter 

angle of yaw (acute angle between the spanwise axis of 
the model and the plane perpendicular to the free-stream 
flow direction) 


specific weight of model material 
Subscripts: 

D based on D 

L based on L 
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1 local value external to boundary layer 

0 at stagnation line of cylinder 

s just behind bow shock 

T at transition 

t stagnation value as determined by isentropic deceleration of 

flow external to boundary layer 


W 

at wall 

on surface of model 

max 

maximum 


A 

based on A 

00 

in free 

stream ahead of how shock 


Superscript: 

' evaluated at reference quantities T' and p' 

A bar over a symbol denotes the average value over the forward half 
of the cylinder unless otherwise noted. 

APPARATUS 

Description of Tunnel 


The investigation was made in a blowdown jet in the Gas Dynamics 
Branch at the Langley Laboratory. The test section is 12 inches wide 
and 13 inches high. The complete system consists of a high-pressure 
settling chamber, a convergent-divergent nozzle, a pressure-tight chamber 
containing the model support system, and a fixed diffuser . The stagnation 
temperature and pressure are controlled automatically at a remote station. 

The air supply is obtained from tanks where the air is stored at a 
maximum pressure of 5 jG 00 pounds per square inch. The air is reduced in 
pressure by two successive regulating stations to the desired stagnation 
pressure, which was between 200 and 500 pounds per square inch for this 
investigation. The specific humidity of the air was less than 1 part 
of water per million parts of air by weight. The air is heated by a 
combination of a Dowtherm heat exchanger and electric heaters. The maxi- 
mum stagnation temperature which can be attained is l,0k0 Fj however, 
for this investigation the highest temperature was limited to about 
100° F and the maximum variation of the stagnation temperature with time 
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was limited to 1° F per minute. The maximum running time available for 
this particular blowdown jet is from 20 to 40 minutes , depending on 
stagnation conditions. 

Mach number distributions, as obtained from total-pressure surveys 
made in a horizontal plane l/4 inch above the nozzle center line, are 
presented in figure 1. This figure shows a slight variation in Mach 
number in the test region; however, no large disturbances appear. The 
average Mach number is 4.15 with a maximum variation of ±0.03 over most 
of the test region. 


Models 

The heat-transfer model was a 24-inch- long cylinder made of electro- 
lytic pure copper. The cross-sectional shape of the cylinder was a semi- 
circle of l/2-inch diameter. The model was supported in the test section 
by a steel bar with the leading edge grooved out to reduce heat conduction 
between the model and its support bar. A sketch of the. model and its 
support is shown in figure 2(b). The model was fastened to and alined 
with the support by means of small screws and pins. The space between the 
model and the edges of the support was sealed with a polysulfide rubber. 

The model was covered and water-cooled by means of a shield (fig. 2) 
during the transient starting conditions . The shield covered the model 
until equilibrium conditions were established in the test section, and 
then it was rapidly retracted by a cable attached to a pneumatic cylinder. 
The removal time was less than 0.01 second. 

The model was provided with end plates (fig. 2) in order to minimize 
aerodynamic end effects . The upstream end plate wan attached to .the model 
and the downstream end plate was attached to the shield. The clearance 
between the model and the shield in the retracted position was filled by 
a contoured block attached to the model. This provided a smooth end plate 
during the time that temperature data were recorded. Slots were cut into 
the model near each end plate and filled with an insulating material to 
minimize heat conduction at the ends of the model. The end plates were 
parallel to the tunnel sidewalls and li inches from them for all the yaw 

angles tested. The aspect ratio of the model (between the end plates) 
thus varied with the yaw angle according to the relation 

1 = 18 
D cos \|r 

The temperature of the model was measured by copper constantan 
thermocouples which were silver-soldered to the rear of the model as 


8 


NACA RM L56EO9 


shown in figure 2(h) . The thermocouples were spaced 1 inch apart along 
the span of the model. Eleven thermocouples were used so that longitu- 
dinal heat conduction and aerodynamic end effects (if any) could he 
evaluated. The model was located in the hlowdown jet so that the center 
thermocouple was at station 0 in the nozzle (see fig. l) for all yaw 

angles except 0°, for which the model was 4^ inches farther down the 

tunnel. 

Two pressure-distribution models were used to obtain data at 0° 
and 45° yaw. These models were complete circular cylinders and had end 
plates similar to those of the heat-transfer model. The aspect ratios 
between the end plates were 4 and 14, respectively. Pressure orifices 
were installed at intervals of 5° over the first 40° and thereafter at 
intervals of 10° and 15°. 


Instrumentation 

The temperature time history of the model was obtained from a 
36-channel recording oscillograph having elements with a sensitivity 
of 12.8 microamperes per inch of deflection. The elements were so posi- 
tioned that ll inches of deflection was obtained. This corresponded to 

a rise of 250° F above cold junction temperature. The galvanometer 
elements were individually calibrated before each run, giving an accuracy 
in absolute temperature of 1° F. However, on the temperature-time 
records for any one channel, the relative temperatures were accurate to 
within l/4° F. 

Tha stagnation temperature was obtained from a self-balancing 
recording potentiometer. Measurements were made with three thermocouples 
installed in the settling chamber. The stagnation temperature is accu- 
rate to within ±2° F. Stagnation pressure was measured on a Bourdon gage 
accurate to ±2 pounds per square inch. The stagnation temperature and 
pressure were also recorded on the oscillograph to coordinate them with 
the temperature time history of the model. 


TEST PROCEDURE AND DATA REDUCTION 
Heat-Transfer and Recovery-Temperature Tests 


The piping system, tunnel, and nozzle were preheated at subsonic 
flow rates, for a period of about 2 hours prior to each test, until the 
inside wall temperatures of the nozzle were up to about 95 percent of 
absolute stagnation temperature. Immediately after the desired preheat 
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temperature was attained, the supersonic run was started. When the tun- 
nel stagnation conditions reached steady values (after 15 to 20 seconds) 
the model shield was retracted and the cooling water was simultaneously 
cut off. The oscillograph records showed that the film of water left on 
the model after retraction of the shield was evaporated in less than 
0.1 second. The end of this time interval was taken as "time zero" for 
the model temperature time history. Sufficiently large water flow rates 
were used so as to insure a uniform temperature distribution along the 
model span at time zero. A typical time history of the model temperatures 
is shown in figure 3 * 

Recovery temperatures were obtained directly from test data by con- 
tinuing the run until the recorded model temperatures reached a constant 
value. This required from 2 to 3 minutes running time, depending on the 
stagnation conditions. 


The equation for the average heat-transfer rate per unit surface 
area to a slender cylinder is 


The temperature derivatives normal to the cylinder axis have been neg- 
lected so that T would be some effective mean temperature at any cross 
section. For this investigation, however, the heat-transfer coefficient 
was evaluated simply by using the measured temperatures for both T 
and in equation (l). Since the thermocouples were located on the 

back of the model, lower indicated heat-transfer coefficients for small 
times would be obtained than if the actual mean temperatures and wall 
temperatures had been measured. By using available solutions for unsteady 
heat transfer to cylinders with a uniform heat-transfer coefficient at the 
surface (ref. 14), the error in heat-transfer coefficients due to thermo- 
couple location was estimated to be between 1 to 2 percent for all heating 
rates encountered in the present investigation, even at the small times 
at which the data were reduced. For the half -cylinder model used in this 

investigation, ^ The values for the thermal properties of copper 

were taken from reference 15, P- 268 . 

The last term on the right in equation (l) represents the net heat 
transfer due to spanwise heat conduction. Since this depends on the 
second derivative of an experimental quantity, it is difficult to evalu- 
ate accurately. Consequently, heat-transfer rates were calculated from 
the data at small times only, when spanwise temperature variations were 


Data Reduction 
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small and the temperature variation with time was large. All heat- 
transfer data presented are for times of less than 4 seconds from time 
zero. For this time interval, the largest value of the conduction term 
for any run was found to he only 3 percent of the total heat input, so 
that, in general, the conduction was neglected. Estimates showed that 
the quantity of heat contributed by radiation was also negligible com- 
pared with the storage term. The time-temperature slopes were read, 
directly from plots such as shown in figure 3(b). 

The mean equilibrium temperature T e was taken as the arithmetic 

average of the measured temperatures for the condition of = 0. This 

at 

occurred after 2 to 3 minutes from time zero, depending on the stagnation 
conditions. Radiation was again found to be negligible since the tunnel 
wall temperatures were close to model equilibrium temperatures because 
of the preheating procedure. 


THEORETICAL CORRELATION EQUATIONS 
Laminar Boundary Layer 


Approximate correlation equations for the heat transfer to yawed 
cylinders can be derived in a general form by the method of reference 8. 
The method is applied here to the forward half of a circular cylinder. 

The average heat-transfer coefficient can be written as 

T" Qtt 

h = — O 

ho Tw - Te 0 


The substitution of the expression for 


V 


T-w - T. 


obtained from equa- 


tion (46) of reference 8 results in a general relation for the average 
Nusselt number, 


hD 


h_ PqAqD 
h Q Mw 


cos 



U Ns/x=0 



1/2 


0.84 



0.4 


( 2 ) 
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The ratio h/h Q for the forward half of a circular cylinder has been 

&T 

evaluated by the method of reference 8 for - 1 with = 0. 

The local flow distributions needed for these calculations were obtained 
from the measured pressure distributions shown in figure ^ and from the 
assumptions that the flow outside the boundary layer is isentropic and 
independent of the spanwise flow. The local equilibrium temperature 
was taken as 


where 


5 ® « Jo + (l - 

T t * 1 V J ®lo 


T l 0 _ 1 + 7 - -g - - 1 K, 2 cos 2 Jf 

TTzZI^J 


(3) 


oo 


The results of these calculations are shown in figure 5 where the ratio 
h/h 0 is plotted against for different values of T w /T^. . The corre- 
sponding theoretical results from reference 8 for = 6.9 and the 
experimental values from reference 6 for = 3-9 are also shown in 

this figure for comparison. According to the method of reference 8 and 
the data of reference 6 , the value of h/h 0 on a circular cylinder is 
practically independent of T V/ /T^_, H*, and \|r, as indicated in figure 5. 


The mean equilibrium temperature T e can be evaluated from equa- 
tion (5 4 ) of reference 8 after noting that 


It 

it 




dx ~ 0.8l 


(5) 


to an accuracy of 2 percent over a large range of M^,. The reason for 
this is that the value of M-^ at any given value of x does not vary 
much with or \|r. (See ref. 8.) Thus, to the degree of approximation 

in the theory, 

^ ■ P + °-8i(i - **) jr 2 


(6) 
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Then from equations (2) and (6) and from figure 5, the heat transfer to 
the front half of a circular cylinder can be easily calculated for a 
large range of stream Mach numbers, yaw angles, and temperatures. 


Turbulent Boundary Layer 

Detailed calculations for the local heat transfer to cylinders with 
turbulent boundary layers are possible by various approximate methods. 
(See, for example, refs. l6 and 17.) These methods involve the use of 
the integral momentum and energy equations together with somewhat 
arbitrary assumptions for the velocity and temperature profiles, skin- 
friction coefficients, and Reynolds analogy. The assumptions are visually 
concerned with the extension of low-speed results on the flat plate to 
higher speeds and/or arbitrary pressure gradients. In some cases, how- 
ever, the skin friction and heat transfer may be evaluated from estab- 
lished equations for low-speed flow without the necessity of extensive 
calculations. Thus, for example, the expression for the skin friction 
on a flat plate in incompressible flow may be extended to supersonic 
flow with heat transfer by evaluating the gas properties at some refer- 
ence temperature T' (ref. 18). The heat-transfer coefficient may then 
be evaluated directly from Reynolds analogy. 

The Reynolds analogy cannot, of course, be used for the calculation 
of local heat-transfer coefficients on a cylinder with a blunt leading 
edge since the local chordwise shear stress is zero at the stagnation 
point. Also, the boundary layer would, presumably, always be laminar 
for at least some small distance downstream of the stagnation point. 
However, if the average heat-transfer coefficient is desired, as in the 
present tests, some form of Reynolds analogy can be used, and it should 
be possible to derive simple correlating equations. 

The expression for the average skin friction over the chordwise 
distance x on a cylinder in incompressible flow is assumed to be of 
the form 

? w , C 1 (7) 

2 1/7 '' 

puc /pu^y 7 

where C-]_ is a constant. For a flat plate, « 0.015 (ref. 19). The 

Reynolds analogy in terms of the average heat-transfer coefficient and 
average shear stress is assumed to be of the form 


— k- - = c 2 

pUooCp 



( 8 ) 
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For a flat plate and Prandtl number of O. 7 O, C 2 = 1.2 (ref. 20). Com- 
bining equations ( 7 ) and ( 8 ) then gives an expression for the average 
Nusselt number, which is 


hx 

k 





( 9 ) 


This result may be extended to compressible flow by evaluating all 
gas properties at some reference temperature T', which is defined here 
as simply 


T ' ■ |( T S + T ») 

The reference density p 1 is also taken as the average value 

P' = |( p s + Pw) 

which can be simplified to 



by noting that = p^ and that p^ p g over the forward section of 

a cylinder. The reference velocity Uoo in equation ( 9 ) would properly 

be some velocity behind the bow wave in supersonic flow. Then, assuming 

that to the first order of approximation the heat transfer is independent 

of the spanwise flow, the reference velocity Uoo would have to be the 

normal velocity component behind the shock Ujj . Equation ( 9 ) then 

s 

becomes 


hx 

k ! 


C l C 2 a \ 


fa'UjrX 


6/7 


or, with - = P- on the forward half of a circular cylinder, the final 

D 4 

result would be written as 


hD 

k' 


C-^c 


(p’u N s D 


6/7 


( 10 ) 
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Hllpert (ref. 15, p. 142) has used an equation of this form for the 
average heat transfer to tubes in a cross flow. 


The average recovery temperature T e for a turbulent boundary layer 
may be calculated by assuming that the local recovery factor r is given 
by the expression 


r = a 1 / 5 = IS 


T_ - T, 


T+ - T n 


Then, for the half-cylinder. 


T + 


4 

si 


irt/4 t 


— dx 


'0 


+ 0.8l(l - 

v 


(ID 


by analogy with equation (6) . 


RESULTS AND DISCUSSION 
Pressure Distributions 


The results of the pres sure -distribution tests are given in figure 4, 
where the "normal" pressure coefficient defined as 



is plotted against 0 . Data were obtained at = I.98 for \|r = 0° 

and at = 4.08 for i|r = 0° and i|r = 44. 7°* The data at = 4.08 
were obtained in the same nozzle described in the section entitled 
"Apparatus" but at a time when the throat spacing of the nozzle was 
slightly different than that used in the heat-transfer tests for which 
the average stream Mach number was 4.15 (see fig. l). The = I.98 

data were obtained in one of the other blowdown jets in the Langley Gas 
Dynamics Branch. Figure 4 shows that the normal pressure coefficient 
increases with increasing normal Mach number (Mjj = cos ijrj over the 

entire surface of the cylinder. 

The dashed lines in figure 4 were computed from a modified 
Newtonian equation for the pressure coefficient (ref. 21), 

Cp = Cp^ cos^ 
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where 



with p, /p taken as the theoretical pressure ratio across a normal 

L ol 00 

shock at the Mach number Mjj . Since the data are in good agreement with 

00 

the Newtonian equation, this equation is used (as in ref. 8 ) for calcu- 


shown in reference 8 were derived from the same data shown in figure 4 . 


The measured heat-transfer coefficients are presented in figure 6 
in the form of a "wall" Nusselt number plotted against a Reynolds number. 
Each experimental point shown represents an average taken over the time 
interval from l /2 second to 4 seconds from time zero for 3 or 4 repre- 
sentative stations along the span of the model. The values for the 
thermal conductivity and viscosity of air were taken from references 15 
and 22, respectively. The wall temperatures for all the data were within 
the range of 90° F to 120° F so that most of the variation in T w /T^- was 
obtained by changing Tj.. The stream Reynolds number P^U^D/p^ was 

varied from 5 X 10 5 to 13 X 10^, corresponding to a maximum variation in 
the wall Reynolds number P 00 U 00 D/p v of 1.8 X 10^ to 4 X 10^. Data were 

obtained at the yaw angles_of 0°, 21.0°, 39-7°, 49*3°, and 58 . 5 °. The 
theoretical variation of hD/k w with p^U^D/Py. for laminar heat trans- 
fer according to equation (2) (with h/h D = O. 67 ) is also plotted in 
figure 6 for these same yaw angles. 

Comparison of the data with this theory shows that at zero yaw angle 
the experimental Nusselt numbers are in good agreement with the theoret- 
ical values, while at all other yaw angles the experimental results are 
considerably higher than the theory and the differences between theory 
and experiment increase as the yaw angle and Reynolds number are increased. 
At zero yaw the effect of T^/T^ on KD/ky., as predicted by equation 2 

is in good agreement with the data except above a Reynolds number of 
3 X 10 5 where the data are 5 percent to 8 percent higher than the theory. 



which 


Heat Transfer and Recovery Temperatures 
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The experimental recovery temperatures are shown in figure 7 as the 
ratio of the measured equilibrium temperature to the stagnation tempera- 
ture plotted against yaw angle- Each point is the average of several 
temperatures measured along the model when the indicated variation of 
temperature with time was negligible. At zero yaw the data are in close 
agreement with the value predicted by equation (6) for a laminar boundary 
layer. As the yaw angle is increased, the data tend to approach the line 
from equation (ll) for a turbulent boundary layer. While no definite 
conclusions can be obtained from figure 7 because of the limited accuracy 
of the data and theory, the trends suggest the possibility that the 
increase in heat-transfer coefficients with increasing yaw angles as 
shown in figure 6 is caused by transition to turbulent flow. 

In order to explore this possibility, all the data are plotted in 
figure 8 in terms of the turbulent correlation parameters Nu' and R' 
as defined in equation (lO). Laminar heat -transfer data from references 5 
and 6, converted to these parameters, are also included in this figure. 
These two sets of data (excluding data from reference 5 at the higher yaw 
angles where end effects may be present), together with the theoretical 
values for laminar heat transfer as computed from equation (2), correlate 
with the data from the present tests for zero yaw. This indicates that 
the heat-transfer data of the present tests at zero yaw are laminar. 

Then, since the data at 21° yaw show the usual trend of transitional 
heat-transfer data, the increases in heat transfer with yaw angle in fig- 
ure 6 are apparently caused by transition to turbulent flow and a forward 
movement of the transition region toward the stagnation line. Also, since 
the data at yaw angles from 39*7° to 58*5° correlate on one line with the 
slope given by equation (lO), the transition has apparently reached its 
maximum forward position at about 40° yaw. That is, in effect, a com- 
pletely developed turbulent boundary layer occurred on the cylinder for 
I4-O 0 to 60° yaw. Note that the product * 0.024, as obtained from 

the experimental results for greater than 40°, is about 90 percent 
higher than the flat-plate value of 0.013- 

Comparison of the present data with Hilpert's data (cited in ref. 15 .> 
p. 142), which is for the average heat transfer to a whole cylinder at 
zero yaw angle in subsonic flow, shows that at corresponding values of R' 
his data are between the laminar and turbulent values shown in figure 8. 
This can apparently be explained on the basis of the variation in local 
heat-transfer coefficient around a cylinder in subsonic flow as measured 
by Schmidt and Wenner (cited in ref. 15, p. l4l). Examination of their 
data shows that, for subcritical Reynolds numbers, the contribution to 
the total heat transfer from the back of the cylinder is about the same 
as the contribution from the front. Hence, for subcritical Reynolds num- 
bers, the average heat-transfer coefficient for the whole cylinder would 
be about the same as the average on the front. As the Reynolds number is 
increased, the heat contributed from the back increases rapidly in propor- 
tion to the contribution from the front since the heat transfer on the 
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front is apparently still laminar. Thus, for supercritical Reynolds 
numbers, the average heat transfer for the whole cylinder is greater 
than the average on the front. Consequently, for supercritical Reynolds 
numbers, the values of Nu' from Hilpert's data might be expected to 
fall between the laminar and turbulent values obtained in the present 
tests. 

The relative magnitude of the increases in heat transfer with yaw 
angle for the present data is shown in figure 9 > where the ratio of the 
average heat-transfer coefficient for the yawed cylinder to the value 
at zero yaw is plotted against yaw angle. This ratio was 

obtained by cross-plotting the data in figure 6 at the wall Reynolds 
numbers p^UojyVw- of 1.8 X 10^ and 4 X 10^, which correspond to the 

stream Reynolds numbers listed in figure 9* Data from references 5 
and 6 at lower Reynolds numbers are included for comparison. The results 
of the present tests show that, at stream Reynolds numbers of about 

15 X lcP, the heat-transfer coefficient at 30 ° yaw is about twice the 
value at 0° yaw. As previously indicated by figure 8, this large increase 
in heat transfer is evidently caused by transition to turbulent flow fol- 
lowed by a forward movement of transition toward the stagnation line as 
the yaw angle is increased. Since the boundary layer is laminar at zero 
yaw, this forward movement of transition with yaw is obviously caused by 
some factor or combination of factors which change with yaw angle, such 
as a length Reynolds number or the three-dimensional nature of the flow 
(that is, the secondary flow) which occurs within the boundary layer as 
the cylinder is yawed. This latter effect has been previously observed 
and studied in detail on sweptback wings and rotating disks in refer- 
ences 12 and 13 and in investigations by Owen and Randall reported in 
papers of limited availability. 

The heat transfer decreases by about 35 percent from the peak value 
for an increase in yaw angle from 35° to 60°, as shown in figure 9* This 
decrease apparently corresponds to the decrease in heat transfer with 
increasing yaw angle that would occur on a cylinder with a completely 
developed turbulent boundary layer, as indicated by the previous discus- 
sion concerning figure 8. 


Effect of Yaw on Transition 

Length Reynolds number .- In general there are perhaps two types of 
length Reynolds numbers which may characterize transition for a given set 
of circumstances. One of these length Reynolds numbers is based on the 
length from the leading edge measured along a streamline and the local 
stream velocity. In the case of a yawed flat plate this "strearawise" 
Reynolds number can be expected to have some significance. (See ref. 10.) 
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As an example, the transition data on a yawed flat plate (fig. 15 of 
ref. ll) is compared with predictions based on this type of Reynolds 
number, which may be defined as 



PlVl L 

Pi 


where p^, V^, and are taken as the local "effective" values out- 

side the boundary layer. These local effective values may be calculated 
by assuming (as in ref. 25) that the flow near the surface, but outside 
the boundary layer, has gone through a normal shock (here computed for 
the component of the stream Mach number normal to the leading edge) 
corresponding to the shock from the leading edge and then reexpanded 
until the local static pressure on the plate is equal to the original 

stream static pressure. Thus, assuming that the value of for 

L 

transition is constant and equal to the value for zero yaw, the stream- 
wise length from the leading edge to transition would be 



following values for the length of laminar run normal to the leading 
edge, Lrp cos v|r : 


deg 

I* 

cos 'ijr 

Computed, 

in. 

Experimental 
(ref. 10), 
in. 

0 

2.72 

2.72 

45 

1.00 

.86 

50 

.80 

.70 

55 

.62 

•57 

60 

M 

.45 


The good agreement between the calculations and the experimental data 
indicates that the concept of streamwise Reynolds number, as defined 
here, can be used to determine the effect of yaw on transition on a 
flat plate. 
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This strearawise length Reynolds number loses much of its significance 
in flows with pressure gradient where the boundary-layer thickness and 
velocity profiles do not depend on the distance from the leading edge 
only. On a yawed blunt cylinder the length L along an external stream- 
line cannot be defined if the starting point is taken as the stagnation 
line ^u 1 = 0) since L then becomes infinite, as can be seen from the 

expression for L, which may be written as 


This indicates that the flow in the neighborhood of the stagnation line 
is largely spanwise, and yet the boundary-layer thickness is constant 
and the velocity profiles are similar in this region. Thus, as is well 
known, a length Reynolds number for transition in this type of flow 
should be based on some local characteristic dimension of the boundary 
layer itself, such as the momentum thickness. 

For incompressible flow on a yawed cylinder a "resultant" momentum 
thickness A may be defined as 


Then, from the method of reference 8, the momentum thickness was calcu- 
lated for incompressible flow on a circular cylinder by using the equation 


for the chordwise velocity distribution. The local resultant velocity 
is then given by the expression 


The results of the calculation are summarized in the following table 




and a momentum-thickness Reynolds number is then defined as 


PV X A 


U 1 = 2%^ sln 9 


V 1 = p U N 2 sin 2 9 + U^sin 2 * 



where the ratio 


is given as a function of i|r and x: 
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*, d eg 

■<#? ~ • 

x = 0.2 

x = 0.4- 

x = 0.6 

X = 0.75 

0 

0.11 

0.22 

0.34 

0.46 . 

20 

.12 

.22 

•33 

•45 

40 

•15 

.22 

.33 

.43 

60 

.21 

.27 

• 34 

.43 


/ pIL-D 

These results indicate that the ratio RW|— — is almost independent 
of yaw angle for x greater than 0.4. Thus, if has any significance 


as a parameter for transition on a cylinder, the value of the stream 
Reynolds number for transition should be about the same for all yaw 
angles. Since this is contrary to the experimental results of Bursnall 
and Loftin (ref. 9), who observed a decrease in the stream Reynolds num- 
ber for transition ( ] from 3.7 X 10^ to 2.1 X ISp for an increase 

\ ^ /T 

in yaw angle from 0° to 60°, the parameter R^ apparently cannot be used 


to characterize the effect of yaw on transition on a cylinder in incom- 
pressible flow. This implies that some mechanism other than the 
Tollmien-Schlichting type of instability is present. Since the corre- 
sponding calculations of R^ for compressible flow have not yet been 

carried out, no definite conclusion as to the effect of viscous insta- 
bility is possible for this case except for small yaw angles where the 
effects of compressibility would still be small. 


Three-dimensional boundary- layer instability .- In an investigation 
by Owen and Randall a parameter X which characterizes the stability of 
the boundary layer on a yawed cylinder or sweptback wing was defined as 


X = 


P&V, 


max 


( 12 ) 


where & is the local boundary- layer thickness and v max is the maximum 

velocity within the boundary layer in the direction perpendicular to the 
local external streamline. The velocity profile in this direction has 
an inflection point and is, therefore, dynamically unstable. An expres- 
sion for X in terms of the stream Reynolds number, yaw angle, and 
chordwise velocity distribution can be derived from an appropriate solu- 
tion of the laminar -boundary- layer equations. For the incompressible 
case where the chordwise flow is independent of the spanwise flow this 
expression is 
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sin \|f 


172 


^cos \|r 


ru 


SM ♦tan 2 * 


U- 


N„ 



(13) 


according to Owen and Randall, where 


U- 


is a calculable function 


N„ 


U n 


of — — . From an analysis of wind-tunnel and flight-test results they 
■*■’00 


showed that, when the value of X reached approximately 125 on swept- 
back wings, characteristic streamwise striations were first observed in 
the agent used for visualizing transition (usually china clay) . These 
striations were interpreted as the traces of small vortices (alined with 
their axes approximately parallel to the external stream) which, according 
to the analysis, are caused by the dynamic instability of the secondary 
flow. When X reached about 175 the flow had become completely turbu- 
lent with transition moving forward to the neighborhood of the leading 
edge. 


In order to compare these results for sweptback wings with the case 
of the yawed cylinder, equation (15) will be evaluated for a circular 
cylinder. The details of the method used by Owen and Randall for calcu- 
lating F are not available; however, this function may be calculated 
from the integral method of reference 8 applied to incompressible flow. 
The result is 



= CA 


3/2 u l 


Ui 


N„ 



(14) 


where C *» 0.045 hy comparison with the approximate formula of Owen 

U-, 

and Randall. Then, the use of — *■— = 2 sin 0 for the circular cylinder 

U Ncc 

in equations (l3) and (l4) gives the expression 



).045 A^ 2 sin 9 cos 0(4 sin% + tan 2 \|r)l (15) 

l/cos 'dr *— 


|fcos ijr 


which is plotted against \|r in figure 10 for 0 = 50°. The maximum 
value of the stability parameter X max occurred at 0 « 5 O 0 for all 

values of* \|r. Using the values shown in figure 10 and the critical 
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Reynolds number for transition from reference 9 for a circular cylinder 
gives values of « 110, 150, and 130 at \|r = 30°, 45°> and 60°, 

respectively. These values Of X T are in the same range as those of 

the Owen and Randall investigation, indicating that transition on the 
yawed circular cylinder in reference 9 is probably initiated by the same 
mechanism as on the sweptback wings . 

Whether these values of X can be used to predict the onset of 
instability or transition in compressible flow is not yet known, although 
data reported in reference 12 for a swept wing at M<x, = 1.6l indicate 
that the instability occurs at a much lower Reynolds number than predicted 
by the Owen-Randall criterion. For purposes of comparison, a rough esti- 
mate of transition Reynolds numbers in the present tests can be obtained 
by extrapolating the data in figure 6 at 21° and 39*7° yaw back to the 
theoretical laminar lines for these yaw angles, and converting the "wall" 
Reynolds number PooUooD/bw to the stream Reynolds number pJJ^D/p^ 

These estimated values of transition Reynolds numbers are then compared 
in the following table with values predicted from equation (15) (or 
fig. 10) with Xij, = 130 (as determined from the data of reference 9 ) • 



^max 

PcoUcoD 

/pjuA 

deg 

JpcPcfi 

K 

from figure 10 

^00 

for X^ = 130 

\ Poo I T 
from figure 6 

21.0 

0.125 

11.0 X 10 5 

3.6 X 10 5 

59-7 

.220 

3-5 x 10 5 

1.8 X 10 5 


The experimental values for the transition Reynolds numbers are seen to 
be about l/3 to l/2 of the values predicted by the stability parameter 
for incompressible flow. In reference 12 the experimental transition 
Reynolds number (not to be confused with the Reynolds number for the 
appearance of vortices) was about l/lO of the predicted value, using 


\nax _ 


125_ 


(11.8 x 10 6 ) 

and using X^ = 175 


1/2 


O.O363 to correspond to their test conditions 


The above comparisons, together with the previous discussion on 
length Reynolds numbers, tend to indicate that the predominating cause 
of transition on a yawed cylinder is the inflectional instability of 
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the boundary layer. If this is the cas e , the transition Reynolds numbers 
obtained in the present tests (with tunnel turbulence and model roughness 
effects present) may reasonably be expected to be of the same order as 
corresponding values for free-f light conditions. Further evidence sup- 
porting this contention is that the inflectional instability of a three- 
dimensional boundary layer is very much more powerful than the viscous 
or Tollmien-Schlichting instability of a two-dimensional boundary layer, 
as shown in reference 13. 


CONCLUSIONS 


Experimental measurements of the average heat transfer on the for- 
ward half of a circular cylinder at a stream Mach number of 4.15 and 

stream Reynolds number range of 5 X 1()5 to 13 X 1(P (based on cylinder 
diameter) resulted in the following general conclusions: 

1. The average heat-transfer coefficients at zero yaw angle were in 
good agreement with values predicted by laminar theory and also correlated 
satisfactorily with other experimental results at lower Reynolds numbers . 

2. At 35° yaw angle the heat-transfer coefficients were 50 percent 
to 100 percent higher than the values at zero yaw. This trend is in the 
opposite direction of that shown by previous data at lower Reynolds num- 
bers where the heat-transfer coefficients at 35° yaw were about 20 per- 
cent lower than the values at zero yaw. The increase in heat transfer 
with increasing yaw angle, as observed in the present investigation, is 
attributed to transition from laminar to turbulent boundary layer which 
apparently occurred primarily as a result of the inflectional instability 
of the three-dimensional boundary layer on the yawed cylinder. 

C* 

3. A further increase in yaw angle from 35° to 60° resulted in a 
decrease in heat-transfer coefficients (from the peak values) of about 
35 percent. These data, when plotted in terms of typical turbulent 
correlation parameters, indicate that the magnitude of this decrease 
corresponds, approximately, to the decrease in heat transfer with 
increasing yaw angle that would be expected on a cylinder with a com- 
pletely developed turbulent boundary layer. 

4. The average recovery temperatures for zero yaw were in good agree- 
ment with the values predicted by laminar theory. The recovery tempera- 
tures at other yaw angles were between an upper limit corresponding to 

a theory for the turbulent boundary layer and a lower limit corresponding 
to a theory for the laminar boundary layer. 
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5. Estimates based on the results of the present investigation 
indicate that the Reynolds numbers (based on free-stream conditions and 
cylinder diameter) for transition on a yawed circular cylinder are about 

14- X 1C)5 and 2 X 10^ for yaw angles of 20° and 40°, respectively, and for 
the Mach number and temperature conditions of the tests. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 27, 1956. 


NACA RM L 56 EO 9 


25 


REFERENCES 


1. Jones, Robert T. : Effects of Sweep-Back on Boundary Layer and Sepa- 

ration. NACA Rep. 884, 1947. (Supersedes NACA TN 1402). 

2. Sears, W. R.: The Boundary Layer of Yawed Cylinders. Jour. Aero. 

Sci., v«l. 15, no. 1, Jan. 1948, pp. 49-52. 

3. Goland, Leonard: A Theoretical Investigation of Heat Transfer in the 

Laminar Flow Regions of Airfoils. Jour. Aero. Sci., vol. 17, 
no. 7 , July 1950, pp- 436-440. 

4. Fluid Motion Panel of the Aeronautical Research Committee and Others: 

Modem Developments in Fluid Dynamics. Vol. II, ch. XIV, sec. 27 O, 

S. Goldstein, ed., The Clarendon Press (Oxford), 1938, P- 631- 

5- Feller, William V.: Investigation of Equilibrium Temperatures and 

Average Laminar Heat-Transfer Coefficients for the Front Half of 
Swept Circular Cylinders at a Mach Number of 6.9- NACA RM L55F08a, 
1955- 

6 . Goodwin, Glen, Creager, Marcus 0., and Winkler, Ernest L.: Investi- 

gation of Local Heat-Transfer and Pressure Drag Characteristics of 
a Yawed Circular Cylinder at Supersonic Speeds. NACA RM A55H31, 
1956. 

7 . Eggers, A. J., Jr., Hansen, C. Frederick, and Cunningham, Bernard E. : 

Theoretical and Experimental Investigation of the Effect of Yaw on 
Heat Transfer to Circular Cylinders in hypersonic Flow. NACA 
RM A55E02, 1955- 

8 . Beckwith, Ivan E. : Theoretical Investigation of Laminar Heat Transfer 

on Yawed Infinite Cylinders in Supersonic Flow and a Comparison With 
Experimental Data. NACA RM L 55 FO 9 , 1955 . 

9. Bursnall, William J., and Loftin, Laurence K., Jr.: Experimental 

Investigation of the Pressure Distribution About a Yawed Circular 
Cylinder in the Critical Reynolds Number Range. NACA TN 2463, 1951. 

10. Bunions , H. W., and Bryson, A. E.: The Laminar-Turbulent Transition 

in a Boundary Layer. Proc. First U. S. Nat. Cong. Appl. Mech. 
(Chicago, 111., 1951), A.S.M.E. , 1952 , pp. 859 - 868 . 

11. Dunning, Robert W., and Ulmann, Edward F.: Effects of Sweep and 

Angle of Attack on Boundary-Layer Transition on Wings at Mach Num- 
ber 4.04. NACA TN 3473, 1955. 


2 6 


NACA RM L 56 EO 9 


12. Scott-Wilson, J. B., and Capps, D. S.: Wind Tunnel Observations of 

Boundary Layer Transition on Two Sweptback Wings at a Mach Number 
of 1.6l. Tech. Note Aero. 2347, British R.A.E., Dec. 1954. 

13. Gregory, N., Stuart, J. T., and Walker, W. S.: On the Stability of 

Three-Dimensional Boundary Layers With Application to the Flow Due 
to a Rotating Disk. Phil. Trans. Roy. Soc. (London), Ser. A, 
vol. 248, no. 943, July 14, 1955, PP- 155-159- 

14. Jakob, Max: Heat Transfer. Vol. 1, John Wiley & Sons, Inc., 1949, 

P- 275- 

15. Eckert, E. R. G. (With appendix by Robert M. Drake, Jr.): Introduc- 

tion to the Transfer of Heat and Mass. First ed., McGraw-Hill 
Book Co . Inc . , 195® • 

16. Squire, H. B.: Heat Transfer. Vol. II of Modem Developments in 

Fluid Dynamics - High Speed Flow, ch. XIV, L. Howarth, ed.. The 
Clarendon Press (Oxford), 1955, PP- 757-855- 

17 . Kalikhman, L. E.: Heat Transmission in the Boundary Layer. NACA 

m 1229 , 1949 . 

18. Sommer, Simon C., and Short, Barbara J.: Free-Flight Measurements 

of Turbulent-Boundary-Layer Skin Friction in the Presence of Severe 
Aerodynamic Heating at Mach Numbers from 2.8 to 7 .O. NACA TN 3391, 
March 1955- 

19. Falkner, V. M.: A New Law for Calculating Drag. The Resistance of 

a Smooth Flat Plate With Turbulent Boundary Layer. Aircraft 
Engineering, vol. XV, no. 169 , Mar. 1943, pp- 65 - 69 . 

20. Colburn, Allan P.: A Method of Correlating Forced Convection Heat 

Transfer Data and a Comparison With Fluid Friction. Trans. Am. 

Inst. Chem. Eng., vol. XXIX, 1933, PP- 174-210. 

21. Penland, Jim A.: Aerodynamic Characteristics of a Circular Cylinder 

at Mach Number 6.86 and Angles of Attack Up to 90°. NACA RM L54A14, 
1954- 

22. Anon. (Compiled by F. C. Morey): The NBS-NACA Tables of Thermal 

Properties of Gases. Table 2.39 - Dry Air (Coefficients of 
Viscosity). Nat. Bur. Standards, Dec. 1950- 

23. Moeckel, W. E.: Some Effects of Bluntness on Boundary-Layer Tran- 

sition and Heat Transfer at Supersonic Speeds. NACA TN 3653 , 1956. 


Stream mach number, 


naca RM L56E09 


27 


Flow direction 



Horizontal distance, inches 


-4 


-2 


-0 


—2 


— 3.5 


Figure 1.- Mach number distribution based on impact pressure measurements 
in the test region at l/k inch above the horizontal center plane. 
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(a) Photograph showing various components. k“92100«l 
Figure 2.- Photograph and sketches of heat-transfer model. 
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support showing shield in place 


(b) Sketches of model. 
Figure 2.- Concluded. 
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(a) Spanwise temperature distribution. 

Figure 3*- Temperature time histories of the model for a typical run. 
i|f = 0°; p Q = 4ll pounds per square inch gage; T Q = 2^9° F. 
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(b) Temperature variation with time at two stations on the model. 

Figure 3.- Concluded. 
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9 , deg 

Figure 4.- ChorcLwise variation of the crossflow pressure coefficient on 
circular cylinders. The theoretical values of Cp^ were computed 

from the normal-shock equations for the corresponding experimental 
conditions. 
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Figure 5.- Comparison of theory and experiment for the ratio of the 
average heat-transfer coefficient on the front half of a yawed 
cylinder to the alue at the stagnation line. 
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Figure 6.- Experimental and theoretical variation of average Nusselt 
number with Reynolds number on the front half of a circular cylinder 
at five different yaw angles. = 4.15 . The temperature ratios 
listed are average values. 



Figure 7»- Experimental and theoretical variation of average recovery- 
to-stagnation temperature ratio on the front half of a circular 
cylinder. Values used in the theory were = 4.15 and a = 0.7- 
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Figure 8.- Comparison of the average heat-transfer coefficients in terms 
of correlation parameters on the front half of a circular cylinder at 
different flow conditions and yaw angles. 
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Figure 9 .- Effect of yaw angle and Reynolds number on the ratio of the 
average heat-transfer coefficient on the front half of a yawed cylin- 
der to the value at zero yaw. 
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Figure 10.- The variation of the Owen-Randall stability parameter with 
yaw angle at 0 = 50° on a circular cylinder in incompressible flow. 
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